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SUMMARY 

The  objective  of  this  study  is  to  develop  a  simple  empirical  estimate 
for  predicting  the  variation  with  depth  of  a  scalar  concentration  field  in 
the  wake  emitted  by  a  submarine  operating  beneath  the  ocean  surface.  This 
study  treats  only  that  portion  of  the  wake  field  for  which  the  initial 
effects  of  sub  generated  turbulence  have  diminished,  and  for  which  the 
effects  of  ocean  turbulence  and  scalar  field  size  are  the  primary  factors 
affecting  the  process  of  scalar  diffusion. 

After  the  wake  is  1  to  2  hours  old,  the  subsequent  wake  diffusion 
is  dominated  by  horizontal  transport  processes,  vertical  diffusion  being 
Inhibited  by  the  effects  of  ocean  stratification.  The  horizontal  diffusion 
coefficient,  K,  has  been  shown  by  Richardson  and  Batchelor  to  be  proportional 
to  the  4/3  power  of  the  diffusion  scale,  t,  K  =  az4^3.  The  constant  of 
proportionality  'a'  is  evaluated  from  near  surface  ocean  diffusion  data  in 
Section  2  of  this  study.  This  parameter  is  also  evaluated  at  a  depth  of 
300  meters  from  the  one  available  piece  of  ocean  diffusion  experimental  data 
in  Section  3.  It  is  shown  that  there  is  an  apparent  decrease  in  ’a'  (or  the 
diffusion  coefficient  K)  of  approximately  2.5  when  one  compares  the  near 
surface  and  intermediate  depth  (300m)  data. 

It  is  shown  in  Section  4  that  this  depth  dependence  of  'a'  can  be 
estimated  by  application  of  the  inertial  subrange  relation  for  the 
diffuslvity,  K  =  Ae1/3*4^3,  where  e  is  the  rate  of  turbulence  dissipation, 
and  A  is  a  universal  constant.  From  the  available  ocean  diffusion  data  and 
from  the  variation  of  e  with  depth,  A  is  shown  in  Section  4  to  be  nearly 
invariant  with  depth.  The  depth  dependence  of  e  indicates  a  factor  of  2.2 
decrease  in  the  diffusion  coefficient,  K,  between  the  surface  and  300  meters. 

A  correlation  of  the  data  of  Webster  for  the  rate  of  dissipation,  e, 
with  depth  shows  that  e  decreases  as  the  depth  to  the  -.7  power.  This 
variation  is  nearly  identical  to  that  of  the  Brunt-Vaisala  frequency,  N, 
below  the  thermocline.  Hence  a  means  is  provided  for  estimating  the  decrease 
of  e  with  depth  in  different  regions  of  the  ocean. 
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1 .  INTRODUCTION 

The  rate  of  diffusion  of  a  passive  scalar,  emitted  by  a  submarine 
into  the  ocean,  is  determined  by  ocean  turbulence  and  the  size  of  the 
scalar  field  once  the  initial  effects  of  sub-generated  turbulence  have 
diminished.  This  initial  sub-dominated  period  corresponds  to  a  diffusion 
time  of  less  than  1  to  2  hours.  As  noted  by  Ko  and  Alber^,  tne  subsequent 
wake  diffusion  is  dominated  by  horizontal  transport  processes,  vertical 
diffusion  being  inhibited  by  the  effects  of  ocean  stratification.  The 
horizontal  diffusion  coefficient,  K,  has  been  shown  by  Richardson  dnd 
Batchelor  (on  a  theoretical  basis)  to  be  proportional  to  the  4/3  power 
of  the  scale  of  diffusion  i 


K  =  at4/3  (la) 

'a'  is  a  constant  with  units  of  [cm2^3]/sec,  with  £  measured  in  cm.  If 
the  scale  of  the  diffusing  substance  falls  within  a  portion  of  the  ocean 
turbulence  spectra  characterized  by  an  inertial  subrange  (with  dissipation 
rate  e),  then  the  dlffusivity  coefficient,  K,  and  its  constant,  'a',  can 
be  related  to  the  dissipation  rate,  e,  by  the  relation 

K  =  Ae1/3*4/3  Ae1/3  =  a  (lb) 

The  parameter,  A,  appearing  in  equation  lb  is  assumed  to  be  a  universal 
constant,  independent  of  e. 

If  It  can  be  shown  that  A  is  indeed  nearly  a  universal  constant,  and 
if  an  Independent  means  Is  available  for  determining  the  variation  of  e  with 
depth  in  the  ocean,  then  it  Is  possible  to  obtain  an  estimate  of  the  variation 
of  the  diffusion  coefficient,  K,  with  depth.  By  estimating  the  depth  dependence 
of  K*  one  can  then  determine  the  change  in  the  concentration  field,  S(t), 
of  a  passive  scalar  emitted  by  a  source  at  various  depths  in  the  ocean.  It 
is  with  the  objective  in  mind  of  predicting  the  variation  of  a  scalar  con¬ 
centration  field  with  ocean  depth,  for  which  this  study  was  initiated. 
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In  Section  2  below,  near  surface  ocean  diffusion  data  (Reference  2) 
are  examined  to  determine  the  best  estimates  of  the  overall  diffusivity 
parameter,  'a'.  This  parameter  is  also  evaluated  at  a  depth  of  300m  from  the 
one  available  piece  of  ocean  diffusion  experimental  data  (Reference  3)  in 
Section  3.  From  these  observations  and  data  on  the  variation  of  the  turbulence 
dissipation  with  depth,  the  universal  constant.  A,  is  determined  in 
Section  4  and  shown  to  be  nearly  independent  of  depth.  The  Conclusions 
of  Section  5  recapitulates  the  important  empirical  and  theoretical 
formulas  necessary  to  calculate  the  variation  of  a  scalar  field  concentration 
with  depth. 
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2.  NEAR  SURFACE  DIFFUSION  COEFFICIENT  *a' 

From  the  data  correlation  of  Okubo^  [Figure  1],  the  'effective' 
variance  0rc2  of  an  instantaneously  released  dye  patch  in  near  surface 
waters,  grows  with  time  to  the  2.34  power,  i.e., 

arc2  =  -Olt2,34  (2) 


where  o-  2 


rc 


-Jf 


re2S(t,re)27rredre 


// 

^  A 


S2irrQdrQ 
e  e 


=  2axay 


r  =  radius  of  circle  with  area  equal  to  that  within  an  isoconcen¬ 
tration  contour  of  value  S 


The  depths  of  the  dye  patches  analyzed  by  Okubo  varied  from  2  to  60m 
beneath  the  surface  of  the  ocean. 


From  simple  diffusion  theory,  the  horizontal  diffusivity  coefficient, 
K,  of  a  circular  two-dimensional  patch,  with  variance  crrc2,  is  given  by 


K  = 


l  dv 
nr 


2 


Hence  the  diffusivity  K  based  on  Okubo 's  correlation  is 


(3) 


K=.0063t1,34  (4) 

To  calculate  K  in  terms  of  a  characteristic  scale,  £,  it  is  assumed  that 

£  =  2/3ax  =  & arc  s  b arc 

Substituting  for  time,  t,  in  terms  of  £ 

arc  =  (2.45)_1£  =  . 104t1 * 1 7 


or 


t  =  (£/.255)1/1,17 


3 


Picture  1.  Variance,  a  2  vs.  Diffusion  Time  Data 
Compiled  byrc0kubo(2) 
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Hence  K  =  .0063t1,34  =  .0063(2/. 255)1,34/1,17  or 


K  =  .032, 


1.15 


(5) 


If  it  is  assumed  that  the  4/3  law  predicts  an  eddy  diffusivity  which 
is  coincident  with  the  Okubo  correlation  for  K  [equation  5]  for  a  dye 
patch  which  is  1  Km  wide,  i.e.,  =  105cm,  then 


*  3,4/3  *  Wo  *  •03,1-‘5 


.03 


(6) 


e  ,18 
match 


=  .00375 


If  *TOtch  =  104cm,  a  =  .0058. 

Hence,  the  effective  constant,  'a',  for  a  4/3  law  diffusivity,  in 
the  depth  range  60m  <  h  <  2m,  is 


(7) 
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3.  INTERMEDIATE  DEPTH  DIFFUSION  COEFFICIENT  'a' 

The  only  available  data  on  turbu‘ ent  transport  in  stable  intermediate 

( 3 )  . 

waters  is  the  dye  diffusion  measurements  of  Schuert'  made  at  a  depth  of 
300  meters  off  Oahu.  Hawaii.  , 


The  4/3  law  diffusivity  constint  ’a  1 300m  cari  be  determ^necl  .from  a  . 

s  kr 

and  can  be  fitted  to  the  theoretical  relation, 


set  of  data  if  the  measured  variance  o:2  is  known  as  a  function  of  time 

rc  i 


•  2  = 


rc 


|)3a3bn3 


obtained  by  integrating 


(8) 


K  =  ab 4/’3orc1*/3  =  l 


(9) 


The  variance  arc2  of  the  dye  pool  at  a  given! time,  t,  c^  be  determined  from  • 
a  semi-log  plot  of  dye  concentration,  S,  vs.  the  square  of  the  effective’ 


radius,  rg,  by  fitting  the  data  to  the  Gauss i pin  relation 

1 

1  r 

S(t,re)  =  S(t,0)exp  '•  \zhTZz) 


(10) 


where  S(t,0)  =  — — 

7TO  v~r 


rc 

M  =  total  mass  per  unit  depth,  of  mixed  layer 

Such  a  plot  is  presented  .in  Figure  2  from  Schuert 's  data  at  a  time 
46.5  hours  after  the  initial  dye  release. +  From  Figure  2  the  variance  a  2 

I  *  ^ 

is  found  to  be 


o  2  =: ,17km2  =  I .7xl09cm2 
rc 

@  t  =  1.68x10s  sec 


(ID 


Unfortunately,  Schuert  was  able  to  obtain  a  complete  concentration  mtap  at 
this  one  time  only.  At  other  periods  the  data  was  insufficient  to  close  the 
isoconcentration  lines. 
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This  data  point  is  plotted  in  Figure  1,  for  comparison  with  the  near 
surface  data  compiled  by  Okubo.  One  notes  that  Schuert's  one  data  point 
falls  below  most  of  the  near  surface  data.  The  value  of  the  constant  'a' 
required  to  fit  a  4/3  law  theoretical  curve  of  variance  vs.  time  (equation  8) 
through  Schuert's  single  data  point,  is 

(12) 


a  ss  1.6x10 

0  300m 
depth 


-3  cm 


2/3 


sec 


4» 

Note  that  the  value  of  a  is  not  to  be  confused  with  the  value  of  a  report 
by  Schuert  of  6.6x10-3.  The  constant  o  appears  in  the  1962  Okubo  theory( 
^  where  K  =  9/8a2r2'3t  and  S  =  (const/t3)exp[-r4' 3/o2t2]. 


» 


8 
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4.  EVALUATION  OF  THE  UNIVERSAL  CONSTANT  A 

In  order  to  calculate  the  assumed  universal  constant  'A'  appearing 
In  the  4/3  law  dlffuslvlty  formula,  (equation  1)  and  to  determine  If  depth 
variations  can  be  adequately  predicted  by  this  relation,  data  on  the 
variation  of  the  rate  of  energy  dissipation,  e,  with  depth  are  required. 

Values  of  the  rate  of  turbulent  dissipation,  e,  have  been  inferred 
from  measurements  of  the  inertial  subrange  spectrum,  E(k),  for  ocean 
currents  using  the  relation 

E(k)  =  (1.4)e2/3k"5/3  (13) 

The  constant  1.4  is  taken  from  isotropic  turbulence  data. 

Measurements  of  ocean  velocity  spectrum  have  been  made  by  Grant, 
Stewart  and  Moillet  (1962)^,  Stewart  and  Grant  (1962)^,  Grant,  Moillet 
and  Vogel  (1963) ^7\  and  Webster  (1969)^.  The  most  recent  and  most 

/  Q  \ 

extensive  measurements  of  ocean  spectra  are  those  of  Webster'  ' ,  who  made 
deep  sea  current  measurements  in  the  Atlantic  ocean  from  moored  buoys  using 
current  meters.  The  depth  of  the  water  at  the  measurement  site  (site  D)  is 
2600  meters. 

A  plot  of  the  measured  dissipation  rate,  e,  (in  units  of  ergs/cm3/sec) 
is  shown  in  Figure  3  as  a  function  of  depth.  The  corresponding  Brunt-Vaisala 
frequency  is  shown  in  Figure  4. 

A  best  fit  of  Webster's  dissipation  data  is  given  by  the  straight 
line  (in  log-log  coordinates) 


cm2/sec3 
=  1  gm/cm3 

Note  that  below  50  -  100  meters,  the  Brunt-Vaisala  frequency 


e  =  4.4xlO~2Z~0,7 
Z  =  depth  (cm)  p 


(14) 
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It  thus  appears  that  the  variation  of  e  with  depth  below  the  thermocline 
is  proportional  to  N 


e  -  N 


(15) 


Using  equation  14,  one  can  estimate  the  near  surface  dissipation 
rate  (at  Z  *  10m)  and  the  intermediate  depth  dissipation  rate,  corresponding 
to  Schuert's  measurements  at  300  meters  (see  Table  I  below).  Using  the 
values  of  the  diffusion  constants,  'a'  for  near  surface  and  intermediate 
depth  diffusion,  the  universal  constant  A  is  then  calculated  by  the  relation 
A  =  a/el/3  (see  Table  I). 


Table  I 

2/3 

Z(m)  e  cm2/sec3  a  -  A 

3.53xl0“4  4x10" 3  .0565 

3.20xl0"5  1.6.10'3  .0505 

K  =  a*4/3  =  A e1/3*4/3 
4  =  ^<YC 

The  interesting  result  found  from  combining  the  data  from  several 
different  sources  is  that  the  value  of  the  universal  constant  A  is  nearly 
constant  with  depth,  within  a  range  of  10  -  15%. 

This  result  is  not  a  direct  proof  of  the  validity  of  the  diffusion 
law 


10m 

300m 


since  the  dissipation  rate  and  diffusion  coefficients  were  not  measured  at 

1 7 ) 

the  same  location.  Further,  the  hot  film  data  of  Grant,  Moil  let  and  Vogel'  , 
taken  from  a  submarine  in  700m  water  indicates  a  much  higher  near  surface 
dissipation  rate 

'Grant,  at.  al  “  10'2  t0  l°'3o»2/^3 
Z  =  15  -  27  meters 
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which  would  yield  a  value  of  A  nearly  3  times  that  given  in  Table  1. 

Whether  this  data  is  a  local  phenomena  related  to  measurement  technique  and 
near  surface  flow  behavior  awaits  further  clarification. 

The  only  data  for  e  extending  from  the  surface  to  at  least  300  meters 
in  depth  is  that  of  Webster.  His  data  shows  a  factor  of  11  decrease  in 
e  with  depth  which  corresponds  to  a  decrease  in  the  diffusivity  factor, 

'a'  (~e 1 ^ 3 )  of  2.22.  The  decrease  of  a  with  depth  based  on  the  Okubo  and 
Schuert  data  is 

asurface/a300m  =  2,5  (16) 
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5.  CONCLUSIONS 

The  effect  of  depth  on  the  apparent  diffusivity,  K,  of  a  passive 
scalar  can  be  ascertained  by  estimating  the  variation  of  the  rate  of 
turbulent  dissipation,  e,  with  depth  and  by  applying  the  relation 


where  l  =  /Co rc 

A  «  .05  -  .055 
a  =  A£1/3 
b  =  /S 

The  data  of  Webster  indicates  for  site  D  in  the  Atlantic,  that 


e  =  4.4xlO"2Z‘°'7cm2/sec3 
Z  =  depth  in  cm 

For  extrapolation  to  othe  areas  of  the  ocean,  a  linear  relation  between 
e  and  N  might  be  used  below  the  thermocline 


eZ  =  N 
eZ=50m  NZ=50m 


The  data  of  Schuert  and  Webster  indicates  a  decrease  in  apparent  diffusivity 
at  300m  depth  of 


300m 


surface 


_  1/3 

*p00m] 

Le10m  J 


«  1/2.2  or  1/2.5 


The  variance  of  a  dye  or  scalar  patch  varies  with  e  and  time  as 


rc 


2  =  ftVAeb4et3  if  K  =  Ae1/3*4/3 
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The  peak  concentration,  S,  of  a  circular  dye  patch  decreases  as 
the  reciprocal  of  the  variance 


S 


max 


--■/D7  gm/cm3 
™rc 


or 

M/D 

A3b1*et3 

For  diffusion  from  a  line  source  or  from  a  continuous  source  moving 
at  a  fixed  speed 

Smax  =  Q**/C(4/3  A)3/2b2e1/2t3/2] 

Thus  the  peak  concentration  for  a  given  time  after  release  will  decrease 

1  /2 

inversely  with  e  for  a  point  source  and  e  for  a  line  source.  Hence, 
peak  concentrations  will  decrease  much  slower  with  time  at  depths  than 
near  the  surface,  and  hence  will  be  significantly  larger  than  near  the 
surface; 


SZ=300m  _  e10m 

-10 

point  source 

SZ=10m 

e300m 

instantaneous  release 

SZ=300m  ..  /e10m 

)  ~3 

line  source  or 
continuous  (sublike) 
source  at  fixed  speed 

SZ=10m 

\E300m. 
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